Ramanujan - Fourier Series and 
the Density of Sophie Germain Primes 



H. Gopalkrishna Gadiyar and R. Padma 
AU-KBC Research Centre, M. I. T. Campus of Anna University, 
Chromepet, Chennai 600 044, India 
E-mail: {gadiyar, padma}@au-kbc.org 

Abstract. A prime p is called Sophie Germain prime if 2j9+ 1 is also prime. 
A formula for the density of such primes is given in a more general setting 
using a new approach. This method uses the Ramanujan - Fourier series 
for a modified van Mangoldt function. The proof remains heuristic as inter- 
change of certain limits has not been justified. Experimental evidence using 
computer calculations is provided for the plausibility of the result. 



Introduction. Ramanujan's work has been a source of inspiration to re- 
searchers in the domain of experimental mathematics [Borwein, Borwein, 
Bailey 1989]. Ramanujan himself used extensive numerical calculations be- 
fore making conjectures and one wonders what he would have done with the 
present technology of computers. Fourier series are an extensively used tool 
across disciplines as a means of extracting useful information from seemingly 
disorganized data. In this paper we wish to bring to the attention of the ex- 
perimental mathematics community the Ramanujan- Fourier series which are 
useful for analyzing arithmetical functions. We use this tool to give a heuris- 
tic "proof (which can be made rigorous if certain limits can be interchanged) 
and give numerical evidence for the formula derived. 

A positive integer p is called a Sophie Germain prime if both p and 2p+l 
are primes, (2,5), (3,7), (5,11), (11,23) for example. The question is: Are 
there infinitely many Sophie Germain primes? The conjecture is: Yes. The 
density of Sophie Germain primes has been conjecured in [Hardy, Littlewood 
1922] following probabilistic reasoning of Brun. Sophie Germain primes are 
of great interest in recent times after the famous AKS algorithm for primality 
testing [Agrawal, Kayal, Saxena, 2002]. If the conjecture about the density 
of Sophie Germain primes is true, then the complexity of the AKS algorithm 
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can be brought down to 0(log n). Sophie Germain primes are the most 
sought after primes for the RSA algorithm as they are robust against the 
Pollard's p — 1 method of factoring [Stinson 1995]. 

In [Gadiyar, Padma 1999], we gave a heuristic proof of the twin prime 
conjecture using the Ramanujan - Fourier series. In this paper, we show 
that the "proof could be extended to a more general class of primes in 
which Sophie Germain prime is a special case. Unlike the earlier probabilistic 
arguments, our approach is analytic. We show that the conjecture regarding 
the density of prime pairs {p,p') satisfying the condition ap' — bp — I, where 
a and b are positive integers, is true upto interchange of certain limits. 

We use the following two tools. 1. The Ramanujan - Fourier series for 
d){Ti] 

Ai(n) — A(n) where A{n) is the von Mangoldt function defined as logp 
n 

if n — p^ 1 where p is prime and k any positive integer and is equal to 
otherwise. 0(n) denotes the number of integers less than or equal to n and 
relatively prime to n. 2. Carmichael's formula for getting the Ramanujan - 
Fourier coefficients for arithmetical functions. 

Numerical evidence is given for various choices of a and b which shows 
remarkable accuracy of the conjecture. 

The paper consists of the following sections. In Section 1, we state the 
main "result", in Section 2, we describe the Ramanujan - Fourier series and 
some properties of the Ramanujan sum, in Section 3, we give the "proof 
of main result and in Section 4, we give the numerical evidence of our main 
result for various choices of a and b. 

1. Main Conjecture. Let a, b and I be positive integers, where (a, 6) — 1. 
Let us ask the question whether there are infinitely many prime pairs (p, p') 

satisfying the equation ap' — bp = 11 In [Hardy, Littlewood, 1922], Hardy 
and Littlewood conjectured the following. Let 'n:{a.b.i){N) denote the number 
of prime pairs (p, p') satisfying the condition ap' — bp = I such that p < N. 



unless {I, a) — 1, {l,b) — 1, and just one of a,b,l is even. But if these 



Then 




(1) 
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conditions are satisfied then 



o log^iV WXp-i ' 

^ p\abl \^ / 

p>2 



where 



(3) 



p>2 

and p denotes a prime. 

Main Result. Let "^(a^iN) = En<iv Ai(n)Ai (^). Then, upto inter- 
change of certain hmits, 



hm -^"ifiabDiN) = < 



'fn.i.Mte), if (a,l) = (b,l) = l 

p>2 ^'^ ^ 

and exactly one of a, b, 1 is even, (•4") 

0, otherwise. 



Note that (1) and (2) follow immediately from (4), see for example [Gadiyar, 
Padma 1999]. 

The twin prime problem corresponds to the case a = b = 1 and / = 2. 
The Sophie Germain prime problem corresponds to the case a — 1, b — 
2 and I — 1. In [Gadiyar, Padma 1999], we showed that the twin prime 
problem is related to autocorrelation and hence to the Wiener - Khintchine 
formula which is used in probability and electrical engineering. We show 
in Section 3 that the heuristic proof that was given in that paper can be 
extended to prove (4). To do this, we need a short description of Ramanujan 
- Fourier series which we give in the next section. 

2. Ramanujan - Fourier Series. Ramanujan in [Ramanujan 1918] showed 
that many important arithmetical functions {a{n)) have an expansion of the 
form 

00 

«W = Z]agCg(n), (5) 

q=l 

where 

n 

2m -n 



C. 



n) = E e^-i^ (6) 



fc=i 

{k,q)-- 
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is called the Ramanujan sum and the a^'s are known as the Ramanujan- 
Fourier coefficients. He obtained such expansions for d{n), <j{n), 4>{n) and so 
on where d{n) denotes the number of divisors of n and a{n) denotes the sum 
of divisors of n. In [Hardy 1921], Hardy proved that the Ramanujan sum is 
a multiplicative function of that is, 

Cgq' = Cq{n)cq\n) If (q, q') = 1 , (7) 
using which he obtained the Ramanujan - Fourier expansion of Ai(n). 

Ai(n) = f:^c,(n), (8) 

where //(gf) is the Mobius function defined as follows: 

( \ _ { (-1)^, li q = PiP2---Pk, p/s are distinct primes , , . 
^^^^ ~ 1 0, otherwise. 

If p is prime, then 

^(")- 't: 

where a\b means a divides b and a/6 means a does not divide b. 

Neither Ramanujan nor Hardy gave a formula for finding the Ramanujan 

- Fourier coefficients which are the back bone of Fourier analysis. This was 
done later by Carmichael [Carmichael 1932]. Let M[f) denote the mean 
value of an arithmetical function /, that is, 

M{f) = lim ^ E /H ■ (11) 

For 1 < k < q, {k,q) = 1, let ek{n) = e^''^^", {n e Af). If a{n) is an 
arithmetical function with expansion (5) then 

a, = ^M(a c,) = ^ i^^l Y.^a{n)c,{n) . (12) 



Also, 



f 1 if k = 

M(e, e^)= ' I (13) 
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3. "Proof of Main Result. Using tlie results in Section 2, we give the 
heuristic proof of (4). Note that A(^^) is defined if and only if a\{bn + /). 
Using the identity 

a— 1 ( ^ •£ \ 

„ . T I 1 IT n \ Tn 

(14) 



2iTi — m 



1 , if a\m , 
, if a /m , 



we write 



lim — *fa&n(-^) 



1 



bn + l\ 1 



n<N 



a-l 

v. 

a la r, 



I 



(k,q)=l {k',q')=l 



lim — y e V" 7a Jaj 

n<N 



(15) 



^ 1-^ ~ ^ ~ ^ /^(g)/i(gO -2..(^+,)i 

^hh h h h mm 

(k,q)=l {fc',q')=i 

(16) 

= S (say), 

where we have used (8), freely interchanged the sums and limits to obtain 
(15), and then used (13) to get (16). We prove that 5* is equal to the R. H. 
S. of (4) as a lemma. 
Lemma. If (a, 6) = 1, then 

' f n.i.«(g) , if (a,l) = (b,l) = l, 



p>2 



and exactly one of a, b, 1 is even 



(17) 



0, otherwise. 

Proof. Since {k', q') = 1 and < j < a — 1, 

k _ fk^^ .\b 
q \q' ) a 
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k' + jq' 

= -> ^ 

q'a 

ki 

— — b, where 1 < ki < q'a and {ki, q') = 1. (18) 
a 

Now, (18) can happen if and only if 

qd2 = q'di , (19) 
for some divisors di and ^2 of o and h respectively and 

k = ^hq 
q'a 

ki^^qd2 



b_ 
d2 



= kif, (20) 

is an integer. Since (a, 6) = 1, this can happen if and only if ^ divides ki. 
Also, from (19), di\qd2 and since (a, b) = 1, So we write q = diqi where 
gi > 1 is an integer. Similarly q' = d2q2 where 52 ^ 1 is an integer. Thus 
from (19), qi = q2. Also /x(q) = n{diqi) 7^ if and only if {di,qi) ~ 1. 
Similarly we have (^2,^2) = 1- 

Let us write ki — ^^2- Since {k, q) — 1 and di\q, we have {di, k) — 1 and 
from (20), k — k2^ and thus (A;2, di) — 1. So we can write 

^fr^tr. i='^ 0(cii)0(ci2)0^(g2) ,fr. ■ ^ ^ 

(dl,'J2) = l (A;2,?2<ild2) = l 

(c(2,g2) = l 

That is, 

(dl.'J2) = l 
(d2,y2) = l 

by the definition of the Ramanujan sum (6). Now by the multiplicative 
property (7) of the Ramanujan sum, since q2,di,d2 are pairwise relatively 
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prime, 

1 ^ fi{di)fi{d2) /i^fe) 

' ^ «SS £ 

(^1,52)=! 
(^2,92) = ! 

(a, 92) = ! 

(6,92)=! 
{ab,q2)=l 

as (a, 6) = 1. Writing the series and sums in (23) as Euler products, we get 
By the property (10) of the Ramanujan sum, if > 1 is an integer, then 

Hence we will assume that (a, = 1 (6, /) = 1 so that 

p|i i^i 
If none of a, b, or / is even, then the product 

So we will assume that one of a,b or I is even. But (a, 6) = 1, (a, Z) = 1 and 
(6, /) = 1 and therefore exactly one of a,b or I is even. We can therefore write 
the infinite product 

p>2 
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Thus S = unless (a, /) = 1, (6, /) = 1 and exactly one of a, 6 or / is even, 
but if these conditions are satisfied, then the value of S as given in (17) is 
got by simplifying (26) using (28). 

4. Experimental Evidence. We give now the compelling numerical evi- 
dence of the main result (4) by varying a and b. We have taken the value of 

C ~ 0.660161816 and Ratio is defined by (R.H.S. of (4))/ ^(°'^j0 W _ 

Example 1. We take a — 1, b — 2, I — 1 which corresponds to Sophie 
Germain primes. In this case, the R. H. S. of (4) = 2C = 1.320323632. 

Table 1 



N 


*(l,2,l)(A^) 


*(l,2,l)(A^) 

TV 


Ratio 


50000 


66130.966133 


1.322619 


0.998264 


100000 


132886.401744 


1.328864 


0.993573 


150000 


200755.416380 


1.338369 


0.986517 


200000 


265612.706085 


1.328064 


0.994172 


250000 


331585.551940 


1.326342 


0.995462 


300000 


394316.641234 


1.314389 


1.004515 


350000 


459668.599011 


1.313339 


1.00531 


400000 


521496.993567 


1.303742 


1.012718 


450000 


588393.432192 


1.307541 


1.009776 


500000 


652614.182933 


1.305228 


1.011565 
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Example 2. We take a = 1, 6 = 10, / = 1. In this case, the R. H. S. of (4) 
= ^ = 1.760431509. 

Table 2 



N 


*(l,10,l)(^) 


^(1,10,1) (A^) 

N 


Ratio 


10000 


17107.791529 


1.710779 


1.029023 


20000 


34210.057148 


1.710503 


1.029189 


30000 


51939.100560 


1.731303 


1.016824 


40000 


70219.348038 


1.755484 


1.002818 


50000 


89934.594398 


1.798692 


0.978729 


60000 


106902.836342 


1.781714 


0.988055 


70000 


123796.944818 


1.768528 


0.995422 


80000 


141470.265879 


1.768378 


0.995506 


90000 


159287.348829 


1.769859 


0.994673 


100000 


177824.093558 


1.778241 


0.989985 



Example 3. We take a = 3, 6 = 5, 1 = 2. In this case, the R. H. S. of (4) 
= = 1.173621006. 

Table 3 



N 


*(3,5,2)(iV) 


^(3,5,2) (iV) 

N 


Ratio 


60000 


69649.061665 


1.160837 


1.011013 


120000 


140371.214304 


1.169770 


1.003292 


180000 


211924.646933 


1.177366 


0.996819 


240000 


282504.323361 


1.177106 


0.997039 


300000 


355072.360724 


1.183578 


0.991587 


360000 


423152.712312 


1.175427 


0.998463 


420000 


496296.973007 


1.181662 


0.993195 


480000 


568659.361599 


1.184709 


0.990640 


540000 


642488.622118 


1.189796 


0.986405 


600000 


712048.221861 


1.186749 


0.988938 
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5. Conclusion. If the step (15) could be proved rigorously which involves 
justification of interchange of certain limits, then a whole class of outstanding 
problems including the twin prime problem and the Sophie Germain prime 
problem could be solved completely. We may say that in a precise sense the 
Ramamijan - Fourier scries for the (refined) von Mangoldt function traps 
the fiuctuations in the distribution of primes. It is hoped that the theory of 
Ramanujan - Fourier series could be developed to study various properties 
of arithmetical functions. Numerical agreement between conjecture and ex- 
periment means that this technique could become a common tool and lead 
to further developments in number theory. 
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